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In this talk I review some of the recent developments in the field of random surfaces and the Dynamical 
Triangulation approach to simplicial quantum gravity. In two dimensions I focus on the c = 1 barrier and 
the fractal dimension of two-dimensional quantum gravity coupled to matter with emphasis on the comparison 
of analytic predictions and numerical simulations. Next is a survey of the current understanding in 3 and 4 
dimensions. This is followed by a discussion of some problems in the statistical mechanics of membranes. Finally 
I conclude with a list of problems for the future. 



1. Outline 

There are many approaches to formulating a 
discrete theory of quantum gravity. In this talk I 
will focus on the dynamical triangulations (DT) 
formulation of simplicial gravity. This was first 
developed in the context of string theory and two- 
dimensional quantum gravity {QG2) where the 
discrete approach has proven extremely powerful 
and actually preceded continuum treatments. By 
now we have considerable confidence in the valid- 
ity of the assumptions underlying the DT model 
in the case of coupling of matter with central 
charge c less than one. This confidence stems 
from the agreement of discrete, continuum and 
numerical results. 

The situation for central charge c greater than 
one was clarified in the last year by David [1]. 
After the Introduction a discussion of this work 
will be the first part of this review. This is fol- 
lowed by a discussion of recent numerical tests of 
David's proposal [2]. 

I will then move on to the fascinating issue of 
the intrinsic fractal geometry of QG2 and QG2 
coupled to conformal matter. In any theory of 
gravity it is natural to ask what effect quantum 
fluctuations have on the structure of spacetime. 
The most basic question we could ask of any 



*Research supported by the Department of Energy U.S. A 
under Contract No. DE-FG02-85ER40237 and by Syra- 
cuse University. 



spacetime with a given topology is: "What is its 
Hausdorff dimension?" Our current knowledge of 
the answer to this question will be reviewed, with 
emphasis on a comparison of analytic predictions 
with recent numerical simulations. 

The success of the DT approach to simplicial 
gravity in two dimensions has inspired several 
groups to tackle large scale simulations of the DT 
discretization of Euclidean Einstein- Hilbert grav- 
ity in 3 and 4 dimensions. Here the situation 
is much cloudier, both theoretically and numeri- 
cally. The current status will be summarized in 
the third part of my talk. 

Any rich idea usually has many unsuspected 
spin-offs. The field of random surfaces and non- 
critical string theory is intimately connected with 
the physics of membranes. Some recent numerical 
results on new phases of anisotropic membranes 
will be discussed as an example of the fascinating 
interdisciplinary nature of the subject of random 
surfaces. 

Finally I present a list of challenging future 
problems that I think are essential for progress 
in the field. I hope that some of them will be 
solved by the time of Lattice 98. 

2. Introduction 

String theories may be viewed as 2-dimensional 
Euclidean quantum field theories with particular 
actions and matter content. From this viewpoint 



they may also be considered as 2d-statistical me- 
chanical models on lattices with dynamical ge- 
ometry. Statistical mechanical models on fixed 
lattices often possess special critical points where 
they are scale invariant. Correlation functions 
of generic matter fields $(x) reflect this scale in- 
variance by transforming very simply under scale 
transformations: ($ (r) $(0)) ~ r~ 2A », where 
A^, is the scaling dimension of $, or equivalcntly 
its anomalous dimension in field theory. Com- 
bining scale invariance with locality leads to the 
much larger symmetry of conformal invariance. 
In two dimensions conformal invariance and uni- 
tarity restrict the possible values of critical ex- 
ponents because the unitarisable representations 
of the associated Virasoro algebra form a discrete 
series parameterized by a single real number - 
the central charge c. The central charge deter- 
mines the effect of scalar curvature on the free 
energy of the model [3] . 

For c < 1 the only allowed values are 
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with m = 2,3,4.. 
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It is also understood, through the classic work 
of KPZ [4], how scaling dimensions of fields are 
modified when the lattice becomes dynamical i.e. 
when the model is coupled to two-dimensional 
gravity. The dressed scaling dimensions A$ are 
determined solely by and the central charge 



c via 
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where the string susceptibility exponent j s de- 
scribes the entropy of random surfaces of fixed 
area A viz: 



Z{A) = e A - A A^~ 3 
and 

1 - h 



(3) 



Js (h) = 2 ^{25 - c+ V(l - c)(25 - c) }(4) 

for a Riemann surface of genus h. Eq.(3) implies 
the singularity structure 
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as A — > A c . The renormalized continuum theory 
is obtained by tuning the cosmological constant 
A to the critical cosmological constant A c . In 
this limit the mean area (A) diverges like l/( A — 
A c ). The linearity of 7 S with genus h given by 
Eq.(4) has the remarkable consequence that the 
partition function including the sum over genus 
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is actually a function not of two variables A and 
H but only of the single scaling combination 



x = (A — A c ) exp 
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For the minimal models of Eq.(l) the scaling vari- 
able is 



(A - A c ) exp 
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The discrete formulation of QG2 dates to 1982 
and has proven to be even more powerful than the 
continuum approach [5] . It is very rare in conven- 
tional field theories for the discrete formulation to 
be more flexible than the continuum treatment 
and the fact that it is so here certainly deserves 
notice. 

In discrete QG2 one replaces the 2d-manifold 
£2) by a set of n discrete nodes {i} and the 
metric g a b(£i,£2) by the adjacency matrix CV,, 
where dj = 1 if i is connected to j and CV, = 
otherwise. The connectivity g, of node i is 



The scalar curvature Ri is 
6 - qi 
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The gravity partition function becomes 
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where Zh, n for pure gravity is the number of dis- 
tinct triangulations (T) of n nodes with genus h 
and Zh, n for models with matter is symbolically 
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for a generic matter field $ coupled to gravity 
with action Sug- The density of states Zh, n is 
best computed by dualizing the triangulation to 
a $ 3 graph and counting the number of distinct 
such connected graphs with n vertices that can 
be drawn, without crossing, on a surface of genus 
h or greater. To incorporate the topology of the 
graph one must promote $ to an N x N matrix 
a la 't Hooft. The full partition function Z(A, jit) 
is, indeed, simply related to the free energy of the 
corresponding matrix model: 
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with the identification 
f 

Jf2 



e 



and g = c 



(15) 



3. The c = 1 Barrier 



As we have seen in the Introduction there is a 
good understanding of how scaling dimensions of 
operators are modified by the coupling of grav- 
ity to conformal field theories characterized by a 
central charge c less than one. In several cases 
the models are exactly solvable and the sum over 
topologies is even possible in the double scaling 
limit N — > oo and A — > A c with the scaling vari- 
able x of Eq. (7) fixed [6] . When the central charge 
c exceeds one the string susceptibility exponent 
7s, according to KPZ [4], becomes complex. This 
unphysical prediction is indicative of an instabil- 
ity in the model. What is the true character of the 
theory for c > 1? The c = 1 point is analogous to 
the upper critical dimension d — 4 in the theory 
of phase transitions and, in fact, there are known 
logarithmic violations of scaling at c=l [7]. 

From numerical simulations for c > 1 the fol- 
lowing picture has emerged. Two basic classes 
of models have been carefully investigated. The 
first consists of bosonic matter fields on dynam- 
ical triangulations. In these models one finds a 
branched polymer BP phase for c large (typi- 
cally c > 5). The BP phase is characterized by 
7s = 1/2 and the proliferation of minimal neck 



baby universes (mimbus). For smaller values of 
the central charge (1 < c < 4) the simulations in- 
dicate an intermediate phase with < 7s < 1/2. 
There is no apparent discontinuity in the vicinity 
of c = 1. The second class of models is multiple 
Ising models on dynamical triangulations. The 
distinct species of spins couple through the dy- 
namical connectivity of the lattice. For n copies 
of Ising model there is a spin-ordering transition 
at a critical temperature. At the critical point 
one recovers a c = n/2 conformal field theory. For 
n > 2, but not too large, 7 S increases smoothly 
from to 1/2 with, again, no discontinuity near 
n = 2 (c = 1). For large n the ordered and disor- 
dered spin phases (both with pure gravity expo- 
nents 7 S = —1/2) are separated by an intermedi- 
ate BP phase with j s = 1/2 and vanishing mag- 
netization. The transition from the magnetized 
pure gravity phase to the disordered BP phase 
is a branching transition with 7 S = 1/2. This 
branching is one sign of the expected instability 
discussed above. Similar phenomena are found 
for the q-state Potts model for q large. A plot of 
7 S as a function of the central charge c from ex- 
isting simulations for two classes of triangulations 
(degenerate 7p and combinatorial 7c) is shown in 
Fig.l. 




Figure 1. 7 S vs c for both degenerate 7x> and 
combinatorial Tq triangulations. 



Could there be a nontrivial infrared fixed point 



governing the critical properties of the c > 1 
models? A renormalization group scheme for ma- 
trix models, roughly analogous to the e-expansion 
about the upper critical dimension, was devel- 
oped sometime ago by Brezin and Zinn- Justin 
(BZ) [8]. A typical example is given by the 2™- 
matrix model formulation of the n-Ising spins 
coupled to gravity. The BZ technique consists 
of integrating out one row and one column of the 
matrix to obtain an TV — 1 dimensional matrix. 
Changing TV is equivalent to changing the string 
coupling constant of Eq.(15). This induces 
a flow in the matrix model coupling which leads 
to definite renormalization group flows. One can 
then look for fixed points. The method is only 
qualitatively correct for the well-understood case 
of c < 1 minimal models but has the advantage 
it can be extended to c > 1. What does it tell us 
about c > 1? 

David's idea was to apply the BZ scheme to ma- 
trix models that include branching interactions. 
These interactions generate microscopic worm- 
holes that connect two macroscopic pieces of a 
Ricmann surface. The action for such models is 
given by 
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The trace-squared (gluing) interaction (with cou- 
pling x) corresponds to branching and is naturally 
generated at second order in perturbation theory 
within a rcnormalization-group framework. Such 
models were first treated by Das et al. [9]. For 
pure gravity it was found that increasing the cou- 
pling x leads to a new critical point (g c , x c ) in the 
(g, a;)-plane separating a large- a; BP phase (7 = 
1/2) from the pure gravity phase (7 = —1/2). 
At the (branching) transition j s = 1/3. This re- 
sult extends to the case of the m — th minimal 
model. In this case the branching critical point 
has 7 S = l/(m+l). Applying the BZ RG method 
to this case David found the RG flows 



0g=9- 6s 2 - 2.ga; 
(3 X = 2x — 3x 2 — 6gx 



(17) 



For c > 1 the only true fixed point is the 
BP fixed point. But the influence of the c < 1 



fixed points is still felt for c in the neighbour- 
hood of 1. For c < 1 the scaling dimension of the 
field x (corresponding to the renormalized cou- 
pling x) is positive and wormholes are irrelevant. 
At c = 1 the branching interaction is marginal 
and the matter-gravity fixed point merges with 
the branching fixed point. For c > 1 wormholes 
are relevant and the matter-gravity and branch- 
ing fixed points become complex conjugate pairs 
off the real axis. As a result there are exponen- 
tially enhanced crossover effects which imply that 
an exponential fine-tuning of couplings 
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Figure 2. Schematic RG flows for n > 2 (c > 1) 
multiplc-Ising models from [1]. 
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is necessary to see the flow to the true BP fixed 
point. Without this fine-tuning flows appear sim- 
ilar to those for c < 1. Fig. 2 shows a schematic 
RG flow for the generalized n > 2 (c > 1) 
multiple-Ising model with branching interactions. 
In the shaded region flows are similar to the case 
c < 1 unless the temperature is fine-tuned. A' 



and A" denote pure gravity fixed points and C 
and C" denote branching critical points. The re- 
gions O and D are ordered and disordered spin 
phases respectively. The line x = corresponds 
to the original discretized model of n-Ising spins 
on a dynamically triangulated lattice. On this 
line there is no spin-ordering transition - only the 
branching transition to the BP phase. 
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Figure 3. 7 S vs p for c = from [2]. 
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Figure 4. 7 S vs p for c = 1 from [2]. 



and p c <~ 0.4 (c = 1). For c = 2 one finds in- 
stead a smooth volume-dependent cross-over to 
the BP phase. Furthermore the results are con- 
sistent with the model being only in the BP phase 
in the infinite volume limit. These interpretations 
are supported by an analysis of the specific heat. 
For c = and c = 1 one sees a definite phase 
transition but for c = 2 no self-consistent crit- 
ical exponents can be extracted from finite-size 
scaling of the specific heat curves. 



This nice picture is consistent with the existing 
numerical results and may also be tested numer- 
ically. A first step was reported at this meeting 
[2] . These authors study the partition function 
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where p is a chemical potential, n m is the num- 
ber of minimal necks and standard mat- 
ter action for multiple Gaussian fields. The cases 
studied are zero (c = 0), one (c = f) and two 
(c = 2) scalar fields. They look for a transition to 
the BP phase at a finite critical chemical poten- 
tial p c - The clearest results come from an anal- 
ysis of 7 S and are shown in Figs. 3-5. For c = 
and c = f 7 S changes sharply from 7 ~ —1/2 to 
7 ~ 1/2 at a definite p c with p c <~ 0.6 (c = 0) 



4. Fractal Dimension 

Our current understanding of the purely space- 
time aspects of QG2 coupled to matter is much 
less complete than our knowledge of the effects 
of gravity on the matter fields and critical be- 
haviour. Of basic interest is the intrinsic Haus- 
dorff dimension of the typical surface appearing 
in the ensemble of random surfaces. The Haus- 
dorff dimension is defined by 



d H = 



log A 
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(20) 



where A is the area of the surface and r is some 
appropriate measure of the geodesic size of the 
surface. There is a considerable variety of alter- 
native definitions of dn- In fact A and r here may 
be replaced by any reparameterization-invariant 
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The fractal dimension for c < 1 models: theory and simulations. 
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Figure 5. j s vs p for c = 2 from [2]. 



observables with dimensions of area and distance 
respectively. 

For the case of pure QG2 d{j is known to be 
4 [10,11]. This result employs a transfer matrix 
formalism to study the evolution of loops on the 
surface. One may also use a string field theory for 
non-critical strings [12] to calculate dn- Perhaps 
the simplest approach is to formulate the theory 
on a disk with boundary length i and to use I as a 
ruler for determining the scaling of both A and r. 
In the string field theory approach an ADM-type 
gauge is chosen in which geodesic distance plays 
the role of time. Together with a choice of the 
string field theory Hamiltonian this determines 
the scaling of geodesic distance r with boundary 
length t to be 



(21) 



for the m — th minimal model coupled to QGi. 
The scaling of A vs. £ may be determined by 



standard matrix model calculations. If we assume 
that the area A scales canonically as I 2 we con- 
clude that 



2m 



A ~ r 

implying = 2m. 
lently written as 
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This result may be equiva- 
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A completely different result is obtained by study- 
ing the diffusion of a fermion with the methods 
of Liouville theory [13]. This gives 



d H = -2— = 2 x 
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where a n is the gravitational dressing of a (n + 
1, n + 1) primary spinless conformal field. 

In Table 1 we list the predictions from these two 
formulae together with the results from numerical 
simulations [14,15]. 

Both the analytic predictions discussed above, 
as well as the exact solution, agree on dn = 4 
for pure QGi. A detailed numerical investiga- 
tion of the fractal dimension, determining both 
the scaling of two-point functions defined in terms 
of geodesic distance and the behavior of the loop 
length distribution function, was presented in this 
meeting [15]. For < c < 1 dn is found to be 
very close to 4, in agreement with earlier simu- 
lations [14]. An example of the excellent scaling 
curves obtained is shown in Fig. 6. The current 
best numerical results therefore do not agree with 
either of the analytic predictions discussed above. 
There is not yet complete agreement, to high pre- 
cision, between the various methods of determin- 
ing dn and so one cannot confidently rule out 
both theories. Further work is highly desirable to 




Figure 6. Scaling fits for the two point function 
for QG 2 from [15]. 



consolidate the result that du for < c < 1 is 
independent of the matter. The subtlety of the 
problem when Ising matter is included may be 
seen in [16]. The scaling of area versus bound- 
ary length depends on the precise order in which 
the infinite-volume limit is taken with respect to 
the tuning to the Ising critical temperature. It is 
possible to obtain dn = 4 if the infinite-volume 
thermodynamic limit is taken with T ^ T c fol- 
lowed by tuning to the spin-ordering transition. 
Finally large-scale numerical results for c = —2 
[17], made possible by recursive sampling of the 
space of graphs for this topological model, are in 
excellent agreement with the Liouville prediction 
dn = 3.56. 

5. Higher Dimensional Simplicial Gravity 

Since the DT approach to quantum gravity is 
very successful in two dimensions it is natural 
to explore its implications in higher dimensions. 
There has, in fact, been considerable effort in this 
direction in recent years. Consider, to begin, the 
case of pure Einstein-Hilbert gravity in D dimen- 
sions for D = 3 or 4. The functional integral to 
be evaluated is 
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where M is the spacetime manifold with topology 
chosen from the class Top and the action S(g) is 
given by 



S(g) = Jd^x^g- (A 
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In the following we will restrict ourselves to 
fixed topology since there is little rigorous un- 
derstanding of the meaning of the sum over of 
topologies in the general higher dimensional case. 
The continuum integral over diffcomorphism- 
incquivalcnt metrics is replaced in the DT ap- 
proach by a discrete sum over all possible cellu- 
lar decompositions (gluings) of D-simplices along 
their D — 1-faces with the simplicial manifold re- 
quirement that the neighbourhood of each vertex 
is a Z?-ball and the DT constraint that all edge 
lengths are fixed. So far the majority of numeri- 
cal results have been limited to the fixed topology 
of the sphere. To be specific let us concentrate on 
the 4-dimensional case of the four-sphere S 4 . 
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Figure 7. 4d Monte Carlo time history of No for 
N 4 = 64,000 and k 2 = 1.278 from [22]. The 
horizontal units are 100 sweeps. 



The free global variables at our disposal are 
the numbers Ni{T) of i-dimensional simplices in 
a given triangulation T (i = 0,1,.., 4). These 
5 parameters are constrained by two Dehn- 
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Figure 8. Correlation of vertex number No with 
singular vertex volume from [23]. 



Sommerville relations 
E ("!)*( 2 l + - 1 i)^( T ) = ° * = 1,2(27) 

i=2fc-l ^ ' 

together with the Euler relation 

4 

X(^ 4 )-E(- 1 ) 1 ^( T ) = 2 - ( 28 ) 
o 

These three relations leave two independent vari- 
ables which we may take to be N 2 and N±. These 
are the discrete analogues of the mean scalar cur- 
vature and the volume. The Einstein-Hilbcrt ac- 
tion Eq.(26) then becomes 

Z[k 2 ,k 4 } = exp[-k 4 N 4 (T) + k 2 N 2 (T)} (29) 

T(S 4 ) 

where is the discrete cosmological constant and 
k 2 is the discrete inverse Newton's constant. In 
practice (almost) fixed N± (volume) systems are 
usually simulated by adding a constraint that re- 
stricts the volume to be near a target volume. 
One is then really approximating the fixed vol- 
ume partition function 

Z{k 2 ,N 4 }= cx P [fc 2 7V 2 (T)] . (30) 

T(S 4 ) 

From extensive numerical simulations it has 
emerged that the system described by the parti- 
tion function of Eq.(30) has two distinct phases. 



For k 2 small (strong coupling) the system is 
crumpled (infinite Hausdorff dimension) and the 
mean scalar curvature (R) is negative. For k 2 
large (weak coupling) the system is elongated 
(branched-polymer like) with Hausdorff dimen- 
sion 2 and positive mean scalar curvature of order 
the volume. In the crumpled phase generic trian- 
gulations contain one singular one-simplex with 
two singular vertices [18,19]. These singular sim- 
plices are connected to an extensive fraction of 
the volume of the simplicial manifold. The local 
volume associated with the singular one-simplex 
grows like V 2 ! 3 while the local volume associ- 
ated with the singular vertices grows like V itself. 
The appearance of singular structures is generic 
to simplicial DT gravity in dimensions D > 4. 
For D > 4 one finds a singular D — 3 simplex 
(with local volume of order V 2 / 3 ) and singular 
sub-simplices (with local volume of order V"). 

It is now clearly established that there is a first 
order phase transition connecting the crumpled 
phase with the branched polymer phase in both 
3 [20] and 4 dimensions [21,22]. This is seen dra- 
matically in the time series of Monte Carlo sweeps 
shown in Fig. 7 [22]. A finite-size scaling analy- 
sis of the variance of the N 2 (or cquivalently iVo) 
fluctuations also reveals a maximum which grows 
linearly with the system volume — a classic signal 
of a first-order transition. It requires both large 
volumes (order 64, 000) and long simulations (or- 
der 10 6 sweeps) to see the first order nature of the 
transition emerge clearly in four dimensions. It is 
also becoming clear that the transition itself is 
closely connected to the formation of the singular 
simplices [23]. This is demonstrated in Fig. 8. 

Many characteristics of the elongated phase are 
analogous to those of branchcd-polymers [24] and 
may be reproduced by simple and elegant statis- 
tical mechanical models of branching [25] . With 
a suitable choice of ensemble this class of mod- 
els also possesses a discontinuous transition [26]. 
Thus the weak coupling phase of the theory may 
reflect more about the combinatorial nature of 
the simplicial DT action than about the nature 
of gravity itself. This remains to be seen. 

The lack of a continuous transition for DT sim- 
plicial gravity in higher dimensions has been a 
deterrent for recent progress in the field. From 



a traditional point of view this absence of a crit- 
ical point implies that the model has no contin- 
uum limit we could call continuum quantum grav- 
ity. At least three viewpoints are possible at this 
stage. It may be that the theory is fundamen- 
tally discrete and that no continuum limit should 
be taken. This viewpoint is advocated in differ- 
ent ways in several other approaches to quan- 
tum gravity such as the causal sets formulation 
of Sorkin [27] and the theory of spin networks 
[28]. Alternatively the theory might only have 
an interpretation as an effective theory valid over 
a limited range of length scales. Finally it may 
be that pure DT simplicial gravity is ill-defined 
but that models with appropriate matter or mod- 
ified measures [29] possess critical points and ad- 
mit a suitable continuum limit. This latter ap- 
proach was discussed in this meeting by Izubuchi 
for QG3 [30] . One modifies the action by adding 
terms that effectively enhance higher curvature 
fluctuations these correspond to changing the 
measure locally by powers of y/g. There is some 
evidence that by tuning this extra term one can 
soften the transition to a continuous one at finite 
volumes. It is very likely though that the effect 
disappears in the infinite- volume limit. This issue 
needs further clarification. 

6. Membranes 

The theory of random surfaces with the addi- 
tion of an extrinsic curvature to the action has a 
direct connection with the statistical mechanics 
of flexible membranes [31]. Physical membranes 
are two-dimensional surfaces fluctuating in three 
embedding dimensions. The simplest examples of 
2-dimensional surfaces are strictly planar systems 
called films or monolayers. These alone are sur- 
prisingly complex systems. But when they also 
have the freedom to fluctuate, or bend, in the 
embedding space it is a considerable challenge to 
determine their physical properties. Two broad 
classes of membranes have been extensively stud- 
ied. Membranes with fixed connectivity (bonds 
are never broken) are known as crystalline or 
tethered membranes. Membranes with dynami- 
cal connectivity (relatively weak bonds which are 
free to break and rejoin) are known as fluid mem- 



branes. 

The simplest action for self-intersecting (non- 
self-avoiding) fluid membranes resembles that 
of the Polyakov string with extrinsic curvature 
(bending energy). Since the beta function for the 
inverse-bending rigidity (where the bending rigid- 
ity is the coupling to the extrinsic curvature) is 
well-known to be asymptotically free at one loop, 
the bending rigidity is irrelevant in the infrared 
and self-intersecting fluid membranes are com- 
monly believed to exist only in a crumpled phase. 

Crystalline membranes, on the other hand, 
have a non-linear coupling between elastic 
(phonon) interactions and bending fluctua- 
tions which drives a phase transition from a 
high-temperature crumpled phase to a low- 
temperature orientationally ordered (flat) phase. 

Much can be learned about membranes by 
applying the techniques and experience of lat- 
tice gravity simulations to these condensed mat- 
ter/biological problems. This is illustrated by 
some recent [32] large-scale Monte Carlo simu- 
lations of anisotropic crystalline membranes - 
these are systems in which the bending or elas- 
tic energies are different in different directions. 
We were able to verify for the first time the pre- 
dicted existence [33] of a novel tubular phase in 
this class of membranes. This phase is interme- 
diate between the flat and crumpled phases - 
it is extended (flat) in one direction but crum- 
pled in the transverse direction. Correspondingly 
there are two phase transitions: the crumpled-to- 
tubular transition and the tubular-to-flat transi- 
tion. Thermalized configurations from the three 
phases of anisotropic crystalline membranes are 
shown in Fig. 9. These simulations employed 
a variety of improved Monte Carlo algorithms 
such as hybrid overrelaxation and unigrid meth- 
ods [34] . An even more challenging problem is to 
incorporate the self-avoidance found in realistic 
membranes [35,36]. 

Finally the full physics of fluid membranes, in 
which dynamical connectivity may be modeled by 
dynamical triangulations, is a problem of limit- 
less challenges which ties together common tech- 
niques in lattice gravity and the burgeoning field 
of soft condensed matter physics [37] . 




Figure 9. The three phases of anisotropic mem- 
branes: (a) tubular (b) crumpled and (c) flat. 



7. Future Challenges 

There are many challenges facing the program 
of simplicial lattice gravity. I will give here a few 
outstanding problems and possible future direc- 
tions. 

• Topology Change 

Most of the numerical simulations in the field 
have been on spacetimes with fixed topology. 
From the functional integral point of view it is 
more natural (and perhaps essential) to allow the 
topology of spacetime to fluctuate. In 2d it is even 
possible to perform analytically the sum over all 
topologies in the double scaling limit. A prelimi- 
nary investigation of a 4c? dynamical triangulation 
model of Euclidean quantum gravity with fluctu- 
ating topology was made some time ago by de 
Bakker [38]. 

• Super symmetry 

Whether or not it has physical relevance it 
is clear that one of the most powerful ideas in 
particle physics/quantum field theory/string the- 
ory at present is supersymmetry. Supersymmetry 
severely constrains the analytic structure of any 
model. Almost no progress has been made in for- 
mulating or simulating supersymmetric simplicial 



gravity or supersymmetric random surfaces. If 
we are ever to make contact with critical or non- 
critical superstrings and recent developments like 
duality relations this will be essential. 

• Classical Limit 

Although simplicial quantum gravity does pro- 
vide a non-perturbative definition of quantum 
gravity in four dimensions there is no understand- 
ing of how classical gravity emerges in the long- 
wavelength limit and indeed of how perturba- 
tive scattering amplitudes are reproduced in this 
framework. 

• Triangulation Class 

At present there is no classification of the pre- 
cise class of graphs that result in KPZ, rather 
than Onsager, critical exponents when matter is 
coupled to a dynamical lattice (QG2). For exam- 
ple it has been shown [39] that MDT (minimal 
dynamical triangulation) models in which the lo- 
cal coordination number is restricted to be 5, 6 
or 7 still result in KPZ exponents. 

• Fractal Dimension 

Is the fractal dimension of all < c < 1 matter 
coupled to QG2 really four? 

• Interpretation of QG4 

What is the correct interpretation of 4c? simpli- 
cial quantum gravity given the first order tran- 
sition from the strong to weak coupling phases? 
What is the proper mathematical framework for 
QG4 and how much analytic progress can be 
made? Recent progress in this direction is re- 
viewed in [40]. 

• Renormalization Group 

Recently renormalization group (RG) methods 
have been developed for pure simplicial gravity 
and for simple matter coupled to gravity [41]. 
There is no rigorous understanding of the prin- 
ciples behind the success of these methods and 
further improvement of the technique is highly 
desirable. Further progress in this direction to en- 



able, for example, the direct computation of the 
beta function for random surfaces with extrinsic 
curvature would be very nice. 

I would like to acknowledge Kostas Anagnos- 
topoulos, Simon Catterall, Marco Falcioni and 
Gudmar Thorleifsson for extensive discussion on 
many issues treated in this talk. 
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